On the interpretation of IceCube cascade events in terms of the Glashow resonance 
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The IceCube experiment (IC) has recently observed 2 cascade events with energies between 1 
and 10 PeV. This energy combined with the fact that no muon-track events are observed may be 
interpreted as a cosmogenic v £ interacting in IC via the Glashow resonance (GR) v e e — > W~ — > 
(hadrons or v e e). We point out a unique, background-free signature of the GR, a single isolated 
muon unaccompanied by any shower activity from the interaction v e e — > W~ — > , and propose 
it as a test of this interpretation. We calculate the event numbers and find that a single such 
event is expected over about a three- year period in IC. We also show that, if event rates remain 
at their current levels then, even with the GR, standard cosmogenic fluxes cannot easily explain 
the observations. Moreover, if muon-tracks remain conspicuous by their absence, then new physics 
needs to be invoked. As example scenarios in conformity with the observations, we calculate event 
rates for neutrino decay and Lorentz-invariance violation. 



INTRODUCTION 

Observation of ultra high energy neutrinos would be an 
important step towards identifying the origin of cosmic 
rays and understanding the physics of nature's most pow- 
erful accelerators. After decades of pioneering work [T], 
the IceCube detector [2] is currently the most sensitive 
tool to study this problem. Interestingly, two high energy 
events, with energy in the range of 1 to 10 PeV, have 
recently been reported, with an expected atmospheric 
background of 0.14 events [3]. They are cascade events, 
implying that they could be caused by a) v e or v T inter- 
acting with nucleons (N) via charged current (CC) in- 
teractions, or, b) v a (a — e, /i, or r) interacting with N 
via neutral current (NC) interactions. A u^N CC inter- 
action would have been accompanied by a long charged 
muon-track, which is not seen. Similarly, a v T N CC in- 
teraction at these energies would have been accompanied 
by a "double-bang" 0]. Since CC interaction cross sec- 
tions at these energies are typically higher than NC cross 
sections by a factor of about 2.5, an interpretation via a) 
is favoured. 

Given the present strong constraints on fluxes below 1 
PeV from IceCube, any future observation in the range 
1 — 10 PeV would likely be due to the Glashow reso- 
nance [5], involving the interaction v e e — > W~ , where 
the W~ decays mainly to electrons and hadrons, thus 
creating cascade events. The required resonant energy 
is E = my / /{2m e ) ~ 6.3 PeV, in the right ballpark vis- 
a-vis the two IceCube events. At energies close to the 
resonance, the cross section v e e — > anything is about 350 
times the standard CC interaction, v^N — > anything, 
while the cross section for P e e — > hadrons is about 240 
times the standard CC value. It is for these reasons that 
the GR is an important factor in interpreting any ob- 
served event in the region 1 — 10 PeV [6]. The GR is also 
characterised by a unique and essentially background- 



free signature (the "pure muon" jBj), which, although 
suppressed, may be observable over a three-year period 
in IC as discussed below. These pure, or contained, and 
isolated muons occur when the resonantly produced W~ 
decays into a muon. Hence only a contained muon-track 
(and nothing else) is seen inside the detector. Unlike v^N 
CC scattering, there is no shower activity, and there- 
fore it is easily distinguishable from "standard" events. 
We will extract the necessary flux of v e to generate two 
events, and apply this number to predict the expected 
number of "pure muons" . Since pure muons are essen- 
tially background-free, even the small number of 0.34 
events per year that we find, indicates that in 10 years 
three or four such events at E = 6.3 PeV can be ex- 
pected and serve both as a straightforward test of the 
GR interpretation of IC's two events, as well as provide 
a unique physical signature of this important but as yet 
undetected Standard Model physics. 

In general, neutrino fluxes from energetic astrophysical 
sources like active galactic nuclei and gamma-ray bursts 
are expected to have a significant v„ + component 
which would manifest itself in the form of long muon- 
tracks in IC. So far, IC has not reported any observa- 
tion of such events at energies around the range spanned 
by the reported cascade events. If indeed, further data 
gathering and analysis confirms that they are absent or 
suppressed in the 1-10 PeV range and above it, then 
this would indicate that the incoming neutrino flavour 
composition has a large v e l(y^ + ratio. A natural 
source of mainly electron-antineutrinos are cosmogenic 
neutrinos [3 jS], which at energies between 1-10 PeV are 
indeed dominated by v e created in neutron decay. This 
interpretation of the two IC events has been put forward 
in Ref. [5]) and we investigate it further in this letter 
both qualitatively and quantitatively. We find that typi- 
cal fluxes in the literature have to be scaled up to conform 
to the observations. As at higher energy the cosmogenic 
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neutrinos contain significant fraction of muon-neutrinos, 
we evaluate the implied event numbers for muon-track 
events using several typical fluxes for cosmogenic neutri- 
nos. These rates lie within the present sensitivity of IC 
(red curve in Fig. [TJ and should be observable in the near 
future as more data gathering and analysis is done. 

If indeed, muon events continue to remain unobserved 
at energies between 10 7 — 10 11 GeV, where they are ex- 
pected to dominate in the case of cosmogenic fluxes, new 
physics may need to be invoked to account for their ab- 
sence. We find that if neutrinos decay on their journey 
to Earth, and if the normal mass ordering is realized, 
then the v e fluxes dominate over the z/„. In contrast, 
for an inverted mass ordering the v e and v e are heav- 
ily suppressed. Additionally, we show that tiny amounts 
of Lorentz violation could again lead to muon neutrinos 
(V M + Dp) being suppressed compared to electron neutri- 
nos (l/ e + V e ). 

In the following section we recapitulate the physics of 
the GR and its consequences for IC. 

ICECUBE AND THE GLASHOW RESONANCE 

We shall assume in this note that the two events at IC 
are generated by the Glashow resonance, v e e — > W~ — > 
(hadrons or v e e). It is straightforward to evaluate the 
flux needed at these energies that leads to two shower 
events due to the resonance at IC. They were observed 
during the IC designated 2011-2012 86 string period 
which had a live time of 353.67 days [3]. Hence, we as- 
sume an exposure period of 1 km 3 -yr in our calculations. 
The event rate is jS] 

Rate(^ e e — > hadrons) = 27r —zN a V c s 

18 

/ I d 
dE v J dy -^-{^e e hadrons) <f>p e (E u ) , (1) 

where N A = 6.022 x 10 23 cm" 3 , V cS 1 km 3 , and the 
factor 10/18 takes into account the number of electrons 
in water per unit volume. The events are integrated over 
the upper half sphere since up-moving D e at these energies 
are attenuated by Earth matter. We find 

E 2 <b Be {E v )\ B =63 PoV = 7.6 x 10~ 9 GeV cm" 2 s" 1 sr -1 . 

(2) 

This value is in a region that will be within the reach of 
IC's detection as it collects more data. 

Should the two events indeed be a result of the GR, it 
will be the first observation a) of the Glashow resonance, 
and 6) of extra-galactic neutrinos. Therefore, it is impor- 
tant to confirm whether the observed events are indeed 
due to this process. One way to check this experimen- 
tally would be to look at the number of "pure muon" and 
"contained tau track" events generated as a result of the 



related processes v e e — > and v e e — > v T r respectively, 
which accompany the GR hadronic interaction and occur 
when the resonant W~ decays to a muon or tau instead 
of hadrons [5] . Such pure muon events would be clearly 
distinguishable from the muon-track events coming from 
the neutrino-nucleon interactions because, originating in 
a leptonic interaction, these would not have the charac- 
teristic accompanying hadronic shower. The cross section 
for this process is related to the one to generate hadrons 
by IV.j.p^/IV^hadrons w 0.156. Thus the event rate 
is small but distinctive in signature. With the flux <I>p c 
obtained in Eq. |2]) we can estimate the pure muon event 
rate at E = 6.3 PeV as 0.34 per year, leading to a sin- 
gle event over a span of about three years. However, the 
"smoking-gun" nature of these signals, and the fact that 
the background (from v^e — > v e fi interactions) is negligi- 
ble [TU], makes them a significant tool for GR detection 
and confirmation in spite of the small rates. 

The next step is to consider sources of high energy 
neutrinos responsible for the GR events. Diffuse fluxes 
with a typical power law will have roughly equal and 
v e content when reaching Earth, and one should expect 
at energies below PeV a sizable number of muon-track 
events, which however, so far, have not been seen. It 
is a useful exercise, therefore, to find neutrino sources 
circumventing the problem. We will show in the next 
section that cosmogenic neutrinos are natural candidates 
for this. 



COSMOGENIC NEUTRINOS 

In this section we shall assume that the origin of the 
electron-antineutrinos is cosmogenic [7J [5], i.e., generated 
by interactions of cosmic ray protons with the cosmic 
microwave photons, creating a A resonance: 

P7 -> A+ -> n + n . (3) 

At the low energies, E < 100 PeV, the resulting flux is 
dominated by the neutron decay, and hence 

$° a :$g M :$g T =1:0:0 and$° Q =0, (4) 

is the initial flavour composition. Here a represents all 
flavours of the neutrino spectrum. Lepton mixing will 
modify this flux decomposition, and at Earth we have 
$„ q = P a p $g , where the neutrino mixing probabil- 
ity is P a (s — Yl \ U a i\ 2 \ Upi\ 2 , with U being the lepton 
mixing matrix. Inserting the current best-fit values of 
the oscillation parameters from Ref. [llj . it turns out 
that the (1:0:0) flux composition is transferred to 
(0.55 : 0.25 : 0.20). Within 3cr ranges of the oscillation 
parameters, the maximal value of the ratio $ Pe /$p is 
4.65. 
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Detailed studies of the spectrum of cosmogenic neu- 
trinos arriving at Earth have been made in the literature, 
and we use the one calculated in Ref. [7 j (ESS) to see if it 
can explain the events seen at IC. We normalise this flux 
so that the v?. flux at ~ 10 PeV produces the two events at 
IC over the appropriate exposure. This requires scaling 
the v e flux, and consequently the fluxes of all flavours, up 
by a factor of about 300 from those obtained in Ref. [7]. 
At higher energies, i.e. 10 8 — 10 10 GeV, the decay chain 
A + — > Tr + n — > e + pev e r> e D fl v fl causes a significant bump 
in the fluxes of the electron and muon neutrinos, and 
muon antineutrinos. As a result, there will be a signifi- 
cant number of shower and muon-track events observed 
at IC at these energies. Thus, with these fluxes, IC should 
see about 15 muon-track standard CC events per year in 
the energy range 10 7 - 10 9 GeV, and about 11 similar 
events over energies 10 9 — 10 11 GeV. 

For comparison, we have carried out similar event rate 
calculations using a different prediction for the cosmo- 
genic neutrino flux [5] . With the FRII flux spectra de- 
scribed therein, we find that IC should have seen 2.2 
shower events, and 3.8 muon-track events per year in the 
energy range 10 5 — 10 7 GeV. 

Evidently, if the CC muon-tracks are not observed in 
the near future, then it will clearly be difficult to rec- 
oncile the observation of the two shower events at PeV 
energies with the lack of characteristic muon-track events 
at higher energies using cosmogenic neutrinos. However, 
if we allow for the effect of non-standard physics during 
the propagation of the neutrinos from source to Earth, 
it might be possible to suppress the muon neutrino flux 
while keeping the v e flux unchanged. If this suppression 
is sufficient, it would nicely reconcile the observation of 
shower events with the lack of muon-track events. In the 
following section we investigate two scenarios which allow 
this: neutrino decay and Lorentz-invariance violation. 



NON-STANDARD PHYSICS 

Neutrino decay. Assuming that neutrinos produced 
at the source decay such that only the lightest mass state 
survives, the flux composition reaching Earth is |12j : 



($ e :$ l ,:<f> T ) = (]U ei \ 2 :\U tli \ 2 :\U Ti \ 2 ) 



(5) 



where Q a = $„ Q + $ Sq and either i = 1 or 3, depending 
on the normal or inverted mass ordering. In case of the 
inverted ordering, since |J7 e 3| 2 — 0.02 <C |£^,t3| 2 — §, 
the muon-neutrino flux will dominate and we will run into 
problems with non-observation of muon-tracks. Hence 
we are led to consider the normal mass ordering |13j . 
With the present best-fit values of the mixing parame- 
ters, one obtains ^e/^ — 5.32. Thus, when the fluxes 
are scaled so that the v e flux at PeV energies explains 
the two shower events seen at IC, the induced suppres- 
sion of the + flux at higher energies leads to about 



4.9 events per year. Varying the neutrino mixing pa- 
rameters over their allowed 3cr range, we find that the 
maximum suppression of the muon neutrino flux with re- 
spect to the electron neutrino flux can be $ e /$ M = 20.96, 
when 6>i2 = 30.59°, 6» 23 = 52.96°, 6 13 = 10.19° with 
Sep = 180°. With the ESS fluxes scaled as explained in 
the previous section, but with the corresponding muon 
fluxes suppressed by this factor, the number of muon- 
track events expected in the energy range 10 7 — 10 11 GeV 
reduces to about 1.3 per year. Assuming these optimistic 
values of the mixing parameters, we see that it is possi- 



ble to obtain a marginally consistent picture (see Fig. la I 



explaining the two shower events at PeV energies via the 
GR mechanism and sufficiently suppressing any muon- 
track events. 

Effect of Lorentz-invariance violation (LV). Viola- 
tions of Lorentz-invariance might show up at low energies 
as a tiny effect in neutrino oscillation |14) . To simplify 
the analysis of LV, we study a model where it affects the 
oscillation of only and v e . The effective Hamiltonian 
in the mass basis is modified to 



H e s(a) 



(6) 



2E, 



which gives an oscillation probability of P [v e 
\ sin 2 26»i2 (l - - f£ cos (2J1L)) , where u - 



: 2 
m 2 — m 1 

IE 



and SI = \fu 



Here a is the LV parameter with di- 



mensions of energy, and we have only retained the leading 
a order correction to the probability. The effect of LV 
with various values of a is shown in Fig. |lb| Clearly, 
even with tiny values of the LV parameter, enough sup- 
pression of the V/j, flux with respect to the v e flux can be 
achieved to evade future IC bounds on the muon-track 
events. For example, with a = 10~ 29 GeV, about 0.18 
muon-track events are expected over energies 10 — 10 11 
GeV per year, which implies that IC would only be able 
to observe a single such event over 5.5 years of expo- 
sure. Such a suppression can, therefore, explain the lack 
of muon-track events at high energies, while still leaving 
the v e flux at PeV energies large enough to explain the 
two observed shower events. 

We thus conclude that with non-standard physics it 
is possible to explain, simultaneously, a) IC's two shower 
events at PeV energies (by means of the GR process in- 
volving the cosmogenic P e ), and b) the non-observation 
of muon-track events at the same detector at higher en- 
ergies. 



CONCLUSIONS 

The Glashow resonance is an important Standard 
Model process which should be considered in any in- 
terpretation of IceCube events seen in the energy range 
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FIG. 1: The effect of non-standard physics, neutrino decay and LV (as a function of the LV parameter a), on the 

and ESS fluxes. The fluxes have all been appropriately scaled so that the v e flux at 6.3 PeV explains the two 
shower events observed at IC. The total Vn + za, flux, when affected only by standard oscillation, is identical to the 
v e flux shown here. For the decay scenario, the range of variation of the and flux as the oscillation parameters 
vary over the 2cr and 3cr range is also shown. The projected sensitivity of IC, in the case of a null-event result during 

the period 2010-2012, is indicated by the red line. 



10 6 — 10 7 GeV. It has a unique and background-free sig- 
nature, the pure muon, which can be a smoking gun for 
this process. Natural candidates for neutrinos at these 
(and higher) energies are cosmogenic neutrinos. We have 
shown that standard fluxes, which have to be scaled 
by an (unexplained) factor of 300, do not explain IC's 
two events, since charged current generated muon-track 
events are not seen. While both the need for scaling 
the flux up and the paucity of muon-tracks may change 
as IC collects and analyses new data, we have pointed 
out that the present observations may, however, be un- 
derstood by invoking the advent of non-standard physics 
during propagation. We have demonstrated this with ex- 
plicit examples of neutrino decay and Lorentz invariance 
violation. 
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